Patterns of Chaos Synchronization 
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Small networks of chaotic units which are coupled by their time-delayed variables, are investi- 
gated. In spite of the time delay, the units can synchronize isochronally, i.e. without time shift. 
Moreover, networks can not only synchronize completely, but can also split into different synchro- 
nized sublattices. These synchronization patterns are stable attractors of the network dynamics. 
Different networks with their associated behaviors and synchronization patterns are presented. In 
particular, we investigate sublattice synchronization, symmetry breaking, spreading chaotic motifs, 
synchronization by restoring symmetry and cooperative pairwise synchronization of a bipartite tree. 



I. INTRODUCTION 

Chaos synchronization is a counter-intuitive phe- 
nomenon. On one hand, a chaotic system is unpre- 
dictable. Two chaotic systems, starting from almost iden- 
tical initial states, end in completely different trajecto- 
ries. On the other hand, two identical chaotic units which 
are coupled to each other can synchronize to a common 
chaotic trajectory. The system is still chaotic, but af- 
ter a transient the two chaotic trajectories are locked to 
each other 0, 0] • This phenomenon has attracted a lot of 
research activities, partly because chaos synchronization 
has the potential to be applied for novel secure commu- 
nication systems 0,01 . In fact, synchronization and bit 
exchange with chaotic semiconductor lasers has recently 
been demonstrated over a distance of 120 km in a public 
fiber-optic communication network In this case, the 
coupling between the chaotic lasers was uni-directional, 
the sender was driving the receiver. For bi-directional 
couplings, when two chaotic units are interacting, addi- 
tional interesting applications have been suggested. In 
this case, a secure communication over a public channel 
may be established. Although the algorithm as well as 
all the parameters are public, it is difficult for an attacker 
to decipher the secret message @, 0, @| • 

Typically, the coupling between chaotic units has a 
time delay due to the transmission of the exchanged sig- 
nal. Nevertheless, chaotic units can synchronize without 
time shift, isochronically, although the delay time may be 
extremely long compared to the time scales of the chaotic 
units. This - again counter- intuitive - phenomenon has 
recently been demonstrated with chaotic semiconductor 
lasers [9, [HI [HI, E2] , and it is discussed in the context of 
corresponding measurements on correlated neural activ- 
ity [U,0[IS[j1. 

Several chaotic units may be coupled to a network with 
delayed interactions. Such a network can synchronize 
completely to a single chaotic trajectory, or it may end 
in a state of several clusters, depending on the top ology of 
the network or the distribution of delay times [17I . 1 18L I lQl . 
[20| . Recently another phenomenon has been reported for 
chaotic networks: Sublattice synchronization. If a small 
network can be decomposed into two sublattices, then 



FIG. 1: Two mutually coupled units. 



the units in each sublattice can synchronize to a common 
chaotic trajectory although they are not directly coupled. 
The coupling of one sublattice is relayed by the chaotic 
trajectory of a different sublattice. The trajectories of 
different sublattices are only weakly correlated, but not 
synchronized [2lj . 

In this paper we want to investigate patterns of chaos 
synchronization for several lattices with uni- and bi- 
directional couplings with time delay. There exists a 
mathematical theory to classify possible solutions of non- 
linear differential or difference equations for a given lat- 
tice [12] • However, this theory does not determine the 
stability of these solutions. But in order to describe phys- 
ical or biological dynamic networks, we are interested in 
stable patterns of chaotic networks. The patterns which 
are discussed in this paper are attractors in phase space, 
any perturbation of the system will relax to these pat- 
terns which move chaotically on some high dimensional 
synchronization manifold. 

Our results are demonstrated for iterated maps, for the 
sake of simplicity and since we can calculate the stabil- 
ity of these networks analytically. But we found these 
patterns for other systems, as well, for example for the 
Lang-Kobayashi rate equations describing semiconductor 
lasers. Hence we think that our results are generic. 



II. TWO INTERACTING UNITS 

We start with the simplest network: Two units with 
delayed couplings and delayed self- feedback, as sketched 
in Fig. [TJ 

For iterated maps, this network is described by the 
following equations: 



a t = (1 - e)f(at-i) + £Kf(a t - T ) + e(l - n)f{b t - T ) 
b t = (1 - e)f(bt-i) + eKf(bt-r) + e(l - k)/(o*_ t ) 



(1) 
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FIG. 2: Phase diagram for a = 3/2 (analytical result). 



where f(x) is some chaotic map, for example the 
Bernoulli shift, 



mod 1 



(2) 



with a > 1. In this case, the system is chaotic for all 
parameters < s < 1 and < k < 1. e measures the 
total strength of the delay terms and K the strength of 
the self-feedback relative to the delayed coupling. 

Obviously, the synchronized chaotic trajectory at = b t 
is a solution of Eq. (HJ . Its stability is determined by r 
conditional Lyapunov exponents which describe the sta- 
bility of perturbations perpendicular to the synchroniza- 
tion manifold. For the Bernoulli map, these Lyapunov 
exponents have been calculated analytically [2l|, [23J , and 
for infinitely long delay, r — » oo, one obtains the phase 
diagram of Fig. [2j 

In region I and II, i.e. for 
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(3) 



the two units are synchronized to an identical chaotic tra- 
jectory a t = b t . Although the two units are coupled with 
a long delay r, they are completely synchronized with- 
out any time shift. For t — > oo, this region is symmetric 
about the line k—\. 

Complete synchronization can be understood by con- 
sidering a single unit driven by some signal St: 



a t = (1 - e)f(a t -x) + ekf{at-r) + H- 



(4) 



If the system is not chaotic, i.e. if its Lyapunov exponents 
are negative, then the trajectory at relaxes to a unique 
trajectory determined by the drive s t . For the Bernoulli 
shift, this region is defined by the inequality 
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< k < — 
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(5) 



and indicated by II + III in Fig. [2 Since the trajectory 
which A relaxes to is uniquely determined by the drive, 



further identical units receiving the same drive would 
relax to the same trajectory. Thus Eq. © defines the 
region where identical units which get the same input 
synchronize to a common trajectory. 
Now let us rewrite Eq. (!]): 



(H = (1 - e)/(ot-i) + e{2n - l)/(a t _ T ) + 

+ e(l - K)f(at-r) + e(l - K)f(b t - T ) , 

b t = (1 - e)/(6*-i) + s(2k - l)f{b t - T ) + 

+ e(l - K)f{bt- T ) + e(l - «)/(a t -r). 



(6) 



Both systems are driven by the identical signal 

S t -r = £(l - k) [/(ot_ r ) + f(b t - T )] ■ (7) 
Comparing (J4j) with (jfj), one finds that for 



hi = 2k — 1 



(8) 



the interacting units © are described by the driven unit 
(j4]). The phase boundary of the driven system, region II 
+ III [Eq. J5J]j and the phase boundary of the interacting 
system, region I + II [Eq. |(3j], are connected with each 
other: With the mapping of Eq. ([8]), one phase boundary 
can be obtained from the other. This is not only true for 
t — > oo but for any value of r. 

Additionally, this mapping does not only hold for the 
Bernoulli shift but for any chaotic system, provided that 
the signal does not change the Lyapunov exponents of 
the driven system. Simulations with other maps, e.g. 
with the Tent map or the Logistic map, confirm the ap- 
plicability of the mapping J8]). Since the slopes of these 
maps are not constant, the Lyapunov exponents depend 
on the trajectories; in order to get trajectories which are 
comparable to the trajectories of two mutually coupled 
units, the driving signal in the driver/receiver setup must 
come from an identical unit c t (which has an increased 
self-coupling due to the lack of external couplings), i.e. 



st-r = e(l - ft)/(c{_ T ) with 
ct = (1 -e)/(c t _i) +e/(c t _ 



(9) 
(10) 



Even for the Lang-Kobayashi equations |24Ll25j|. which 
describe coupled semiconductor lasers, we found a good 
agreement for the mapping JH]), see section IVIll 

Now let us come back to the iterated equations |T]). 
Let us assume that we record the synchronized trajec- 
tory at = bt of two interacting chaotic units, regions I 
and II of Fig. [2j Now we insert the recorded trajectory 
bt into Eq. (Q]). How will at respond to this drive? We 
find that in region II the unit A will synchronize com- 
pletely to the recorded trajectory bt, whereas in region I 
the unit A does not synchronize. Although the two inter- 
acting units A and B do synchronize, the unit A does not 
follow the recorded trajectory in region I. This shows that 
bi-directional interaction is different from uni-directional 
drive. 
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FIG. 3: Ring of 4 units. 



III. SUBLATTICE SYNCHRONIZATION 




FIG. 4: Sublattice synchronization in a triangular lattice with 
periodic boundaries. The double lines signify bi-directional 
couplings. The self-feedback is not drawn to simplify the il- 
lustration. 



The response of a single chaotic unit to an external 
drive, Fig. [21 determines also the phase diagram of a 
ring of four chaotic units. Additionally, it shows a new 
phenomenon: sublattice synchronization |2ll |. Consider 
the ring of four identical units of Fig. [3j The dynamics 
of unit A is described by 



a t = (1 - e)/(a t _i) + enf{a t ~ T ) + 
+ e(l - k) 



\fih-r) + \f{dt-r) 



(11) 



the dynamics of B, C and D are defined analogously. Note 
the additional weighting \ for the external coupling (to 2 
neighbors) which causes the total strength of the external 
coupling to remain e(l — k). 

Obviously, the two units A and C receive identical in- 
put from the units B and D. Consequently, they will re- 
spond with an identical trajectory in the regions II and 
III of Fig. [21 The same argument holds for the two units 
B and D, which both get the same input from A and 
C. That leads to sublattice synchronization in region III 
of Fig. [21 A and C have an identical chaotic trajectory 
and B and D have a different one. Although there is a 
delay of arbitrary long time of the transmitted signal, 
synchronization is complete, without any time shift. The 
synchronization of A and C is mediated by the chaotic 
trajectory of B and D and vice versa. But the two tra- 
jectories have only weak correlations, they are not syn- 
chronized. Numerical simulations of the Bernoulli system 
with small values of r have shown that there is no gen- 
eralized synchronization either [2ll . l26l . [27J • 

Sublattice synchronization has been found for rings 
with an even number of units, for chains and also for 
other lattices which can be decomposed into identical 
sublattices [2l|. For example, the lattice of Fig. [4] can be 
decomposed into three sublattices. For some parameters 
of the Bernoulli system we find sublattice synchronization 
with three chaotic trajectories. Again, the synchronized 
units are not directly coupled, but they are indirectly 
connected via the trajectories of the other sublattices. 

The sublattice trajectories, described before, are sta- 
ble, i.e. the conditional Lyapunov exponents which de- 
scribe perturbations perpendicular to the synchroniza- 
tion manifold are negative. Even when the system starts 



from random initial states it relaxes to the state of sublat- 
tice synchronization (region III of Fig. [21 for the ring of 4 
units) . The chaotic trajectories of sublattice synchroniza- 
tion may be depicted as g ? . Note that this structure 
does not break the symmetry of the ring: The statistical 
properties of the chaotic trajectory of A is identical to 
those of B. 

However, there are other solutions of the dynamic 
equations, as well. These solutions are classified accord- 
ing to the theory of Golubitsky et al [22]. For example, 
for the ring, the state g g is a solution as well. Such 
a state breaks the symmetry of the lattice. But we find 
that this state is unstable. Any tiny perturbation relaxes 

to the states g a m re gi° n HI; a a m re gi° n H an d d c 
outside of II and III. In fact, we have never found a sta- 
ble state which breaks the symmetry of the lattice. Only 
when the couplings in the two directions of the ring are 
different, the states g g and ^ | are stable in some (dif- 
ferent) parts of the parameter space, see Fig. [H Hence 
we postulate that spontaneous symmetry breaking is not 
possible for finite lattices of chaotic units. 



IV. SPREADING CHAOTIC MOTIFS 

The response of a chaotic unit to an external drive, 
Fig. [21 points to another interesting phenomenon. Con- 
sider a triangle of chaotic units with bi-directional cou- 

Choose the parameters 



plings as sketched in Fig. 6(a) 



such that the triangle is completely disordered, but each 
unit has negative Lyapunov exponents when it is sepa- 
rated from the two other s. (Both conditions are fulfilled 
in region III of Fig. 6(b) , which shows analytical results 
for the Bernoulli system.) When we record the three time 
series at, bt and c* we find three different weakly corre- 
lated chaotic trajectories. Now we feed the two trajecto- 
ries bt and Ct into an infinitely large lattice of identical 



units with uni-directional couplings as shown in Fig. 6(a) 



Each unit receives two input signals from two other units. 
But since all Lyapunov exponents are negative, the sys- 
tem responds with the three chaotic trajectories a t , bt 
and c t . Although the units of the initial triangle are 
not synchronized, their pattern of chaos is transmitted 
to the infinite lattice. All units of the same sublattice 
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(a) Two different coupling 
parameters rei and K2- 
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(b) Phase diagram for Bernoulli system with 
a = 3/2 and e = 11/20 (analytical result). The 
regions of no synchronisation are labeled with a 
dash (— ). 

FIG. 5: Ring of 4 units with two different coupling parame- 
ters. 



are completely synchronized without time shift, although 
the coupling has a long delay time r. Hence the chaotic 
motif, three weakly correlated chaotic trajectories, can 
be imposed on an arbitrarily large lattice. Note that the 
time for spreading a motif on a large lattice increases 
only linearly with the number of units because of the 
unidirectional couplings. 

For some p arame ters n and s, namely in regions I 
and II of Fig. |6(b)[ the three units of the triangle are 
completely synchronized. In region I, only the three 
units are synchronized while the other units remain un- 
synchronized. In region II, the other units, too, get syn- 
chronized to the triangle (because all Lyapunov expo- 
nents are negative), so the whole lattice is completely 
synchronized. 

The phenomenon of spreading motifs is not restricted 
to a triangle. For example, a ring of 6 mutually cou- 
pled units has also a region in the corresponding (e, n) 
phase diagram where there is no synchronization, but 
where isolated units have negat ive Lyapu nov exponents, 
comparable to region III in Fig. |6(b) [21]. Consequently, 
the pattern of six chaotic trajectories can lead to sub- 
lattice synchronization with 6 different sublattices on a 
corresponding lattice of chaotic units with short-range 




(a) Synchronization pattern (Sublattice 
synchronization) of region III. 
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(b) Phase diagram for Bernoulli system, 
a = 3/2 . Analytical result combining the 
analytical synchronization region for a ring of 3 
units ( — ) and the region where identical units 
which are driven by an identical signal, 
synchronize (- -). 

FIG. 6: Triangle (three bi-directionally coupled units) with 
a uni-directionally attached infinitely large lattice. Double 
lines signify bi-directional couplings whereas arrows show uni- 
directional couplings. The self-feedback is not drawn to sim- 
plify the illustration. 



unidirectional couplings. 



V. SYNCHRONIZATION BY RESTORING 
SYMMETRY 

In general we expect that the larger the network is, 
the smaller the region in the parameter space is where 
the network synchronizes. For example, a ring of N = 6 
units has a smaller region of synchronization than the 
region II and III of Fig. [2] for N = 4. With increasing 
N (and the slope a being held constant) synchronization 
finally disappears completely jll|. However, we found 
a counterexample where adding a unit restores synchro- 
nization. 

" The 



Consider the chain of 5 units shown in Fig. 7(a 
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(a) Chain of five units with two different time delays 
and without self-feedback. 

(b) The last unit has been removed. 

FIG. 7: Symmetric and asymmetric chain without self- 
feedback. 



coupling to the two outer units has a longer delay time 
than the internal couplings. The re is no self- feedback, 
k = 0. Now remove unit E, Fig. |7(b)[ and rescale the 
coupling to unit D. (Unit D doubles the input from unit 
C to compensate for the missing second neighbor.) In 
this synchronized solution does not exist. Nu- 

merical simulations of the Bernoulli system and the laser 
equations show high correlations between units A and C 
with time shift A = n — r 2 , and between B and D with 
zero time shift, but the correlation coefficient does not 
achieve the value one. On the other side, if we add unit 
E we restore the symmetry of the chain. In this case we 
find sublattice synchronization with time shift between 
the outer units and the central one: 



a t = e t = c t _ A ; 



bt = d. 



(12) 



If n is greater than 75, the central unit is earlier than 
the chaotic trajectory of the outer ones, it leads, whereas 
for the opposite case it lags behind. 



VI. COOPERATIVE PAIRWISE 
SYNCHRONIZATION 

Is it possible to synchronize two sets of chaotic units 
with a single coupling channel? In fact, we found such ex- 
amples where two sets of chaotic units are bi-directionally 
connected by the sum of their units, as indicated in Fig.[8l 
There are 2N units, i.e. the number of units on each side 
is N. Each side is the mirror image of the other. Only 
one single bi-directional signal composed of all N signals 
from each side is driving the other side, and this leads to 
cooperative pairwise syn chron ization. 

In the first setup, Fig. 8(a) all units are identical, but 
the delay times of their couplings are different. The units 
have pairwise identical delay times, i.e. A& and B& have 
a coupling delay time 2-Tfc + r which is enforced by a 
self- feedback with delay time t s = 2tu + r. Hence, for 
one pair, N = 1, we obtain the phase diagram of Fig. 
where the two units are completely synchronized in re- 
gions I and II. For a large number N of units, each A 
unit receives the identical signal 



Si 



1 N 




(a) There are N different delay times which 
are pairwise identical for one unit of side A 
and one unit of side B. 




(b) There are N different shifts in the 

Bernoulli map which are pairwise 
identical for one unit of side A and one 
unit of side B. 

FIG. 8: Each unit on one side is coupled to all units of the 
other side. 



and vice versa. Note that the unit A& receives only a 
weak signal of the order 1/N from its counterpart B&. 
Nevertheless, we find that the network synchronizes to 
a state of pairwise identical chaotic trajectories, a,k,t — 
bk,t',k = 1, . . . , N. For the Bernoulli system, the region of 
pairwise synchronization is similar to region II of Fig. [2 
There is no synchronization among units of the same side. 
Each unit receives the sum of all chaotic trajectories, but 
it responds only to the tiny part which belongs to its 
counterpart. The synchronization is a cooperative effect. 
As soon as a single unit is detuned, the whole network 

loses synchronization. 

The second setup, Fig. |8(b)| is similar to the first 
one but allows analytical calculations. In this second 
setup, all delay times are identical, while each unit of 
one side has a different shift (3k, k — 1,...,N, in its 
shifted Bernoulli map, see Eq. (fl4|) : the shifts are pair- 
wise identical for unit Afc and the corresponding unit Bfc 
so one side is again the mirror image of the other side. 
The shifted Bernoulli map is defined by 



fp(x) = [a(x + /?)] mod 1 
Each A unit receives the signal 



1 N 

s t = e(l - k)—22 f 0k (b k 



-r) 



(14) 



(15) 



fe=i 



(2T k +r)) 



(13) 



fc=i 



and vice versa. Analytical calculations are possible for 
this second setup because the delay times are equal - in 
contrast to the first setup - and the different shifts in 
the map do not hamper the calculations. If the Bernoulli 
maps were not shifted, the units would all be identical. 
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Then in regions II and III of Fig. [2j all units of the 
same side would synchronize, aij = fl2,i = • • • = CLN,t\ 
b\,t = b 2 ,t = ■ ■ ■ = bN,t because they receive the same 
input. Then one would effectively get two coupled units 
A and B (which synchronize in regions I and II), so in 
region II all 2N units would be synchronized. For shifted 
Bernoulli maps, the units of the same side are not allowed 
to synchronize by definition. Nevertheless, the stability 
analysis regarding the linearized equations is not affected 
by the different shifts, meaning that the same perturba- 
tions still are damped in the same regions; the only dif- 
ference is that due to the shifts, most of the trajectories 
are not allowed to come close together. Only the pairs of 
corresponding units can synchronize, au,t — &fe,t, and so 
they do in region II. 



VII. ANALOGY OF LASER EQUATIONS TO 
ITERATED MAPS 



Besides iterated maps, we considered the Lang- 
Kobayashi equations, which describe the dynamics of 
semiconductor lasers optically coupled to their own 
or/and to the light of other semiconductor lasers. We 
used them for simulations in the following form, accord- 
ing to Ref. [H: 

-Eoj(t) - -G Nnj {t)E^(t) + ^-^ 

+ A Eo,j (t - t) cos[lu t + 4>j (i) - <f>j (t - t)] 

^lasers 

+ G E u 'j-k E o,k(t - t) cos[cj t + 4>j(t) - <pk(t- t)] , 



k=i 



(16) 



x E 0tj (t-T) ^ 



JV,, 



Wj,k- 



k=l 



E 0}k {t-r) 
EoAt) 



sinjwoT + <j>j(t) - <f>j(t - t)] 
sin[w T + (j)j(t) - cj> k (t - t)) , 



(17) 



-n 3 (t) = (p-l) 7 iV sol -7n i (i)-[r+G Ar n J -(i)] J B 2 J (i) , 

(18) 

where Eoj(t) and (f)j(t) are the amplitude and the 
slowly varying phase of the electric field Ej(t) = 
Ea,j(t) exp{i[woT + <t>j(t)]} and rij(t) is the carrier num- 
ber above the value for a solitary laser, j = 1, . . . , N\ aseis . 
The strength of the self-feedback is determined by A, 
while the strength of the external coupling is defined 
by a. The equations above cover the general case 



of Masers coupled semiconductor lasers, where the net- 
work structure and coupling strengths are determined by 
the weightings Wj,fe, which are normalized so that the 
strength of the total input for each laser is the same, 
12k=r™ w j, k = 1 f° r each laser j (28[. The parameters 
•^soi) Grp_r, aief, 7, r, p, uiq and C sp are chosen according 
to Ref. (H. 

In a very simplified form, the equations above read 



dt' 



— [internal dynamics] + A [self- feedback] 



cr[external coupling] . (19) 



Comparing equations HJ and lfl9l) yields a relation be- 
tween the parameter space of the Lang-Kobayashi equa- 
tions, {A, a}, and the parameter space of the maps, 
{e, k}, if the two following conditions are considered: 

(i) The ratio of the self-feedback to the external cou- 
pling should be the same in both cases. 

(ii) The ratio of the time-delayed terms to the internal 
dynamics should be the same in both cases. 

These two conditions yield the following transformation: 



A 



A + a 



A + CT 



(20) 



The second condition (ii) is not properly defined because 
the Lang-Kobayashi equations lfl9|) - in contrast to the 
iterated equations (U) - are differential equations and the 
internal dynamics is not only the first term on the right- 
hand side of Eq. ([I9| but is also contained in the current 
state x. Therefore the denominator v in (J20J) is not given 
and has to be chosen reasonably. We took v = 180 ns" 1 
so e is between and 1 [because the sum A + a (the 
strength of the re-injected light) should not exceed the 
value of v = 180 ns _1 in our case]. 

In order to measure synchronization, we averaged the 
amplitudes of the electric fields, Eqj, in 1 ns intervalls 
and calculated the un-shifted, A = 0, cross correlation 
function defined by 



C 



A 



(x t ■ yt-A) - (x)(y) 



(21) 



A cross correlation over 0.99 can be regarded as synchro- 
nization in the case of lasers. 

Figure [9] shows both the parameter transformation, 
Eq. (j20l) . from {A, a} to {e, k} and the mapping, Eq. ((HI), 
from two mutually coupled units to driven ones. From 
the comparison of this figure with Fig. El the analogy of 
laser equations to iterated maps can be seen. Besides the 
analogy of features described in section II, one can also 
find this analogy for the features described in sections 
EH E] and ED 
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PIG. 9: Phase diagrams for the Lang-Kobayashi laser equations. Every circle/point represents one simulation. Open circles 
(o) show a cross correlation C > 0.99, which can be regarded as synchronization in the case of the laser equations. Figure (a) 
shows results for two mutually coupled lasers before applying transformation ll20l . wheras (b) shows the same diagram after 
the parameter transformation. The synchronization region of two mutually coupled units (b) can be mapped [using Eq. ([8])] to 
a region (c) which is similar to the synchronization region of a driver /receiver system (d). After the parameter transformation 
lf20]l . the synchronization regions for the Lang-Kobayashi equations, figures (b), (c) and (d), look very similar to the ones for 
the Bernoulli maps, Fig. [2] 



VIII. SUMMARY 

Small networks of chaotic units with time-delayed 
couplings show interesting patterns of chaos synchrony. 
These patterns are stable attractors of the network dy- 
namics. 

The phase diagram of two and four units with mu- 
tual couplings has been related to the properties of a 
single driven chaotic unit. Two interacting units with 
self-feedback can synchronize completely, without time 
shift, even if the delay time is extremely large. When the 
chaotic trajectory of two interacting units is recorded and 
used to drive a single identical unit, it turns out that the 
driven unit synchonizes only in a small part of the phase 
diagram. Hence interaction is different from drive. 

Sublattice synchronization is found for lattices which 



can be decomposed into a few sublattices. Each sublat- 
tice is completely synchronized, but different sublattices 
are only weakly correlated. Synchronization is relayed 
by different chaotic trajectories. The trajectories of each 
sublattice have identical statistical properties. Thus the 
symmetry of the lattice is not broken. There are solu- 
tions of the dynamic equations which break the symme- 
try of the lattice. However, we always found that these 
solutions are unstable. Hence we postulate that stable 
patterns of chaos synchrony possess the symmetry of the 
corresponding lattice. 

Synchronization depends on the symmetry of the net- 
work. When the symmetry of an asymmetric chain is 
restored by adding units, sublattice or complete synchro- 
nization is restored, too. 

Finally, a bi-partite network, where the two parts are 
coupled by a single mutual signal, shows pairwise com- 
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plete synchronization, whereas the units of each part do 
not synchronize. Each unit responds to the weak con- 
tribution of its partner in the other part of the network. 
Pairwise synchronization is a cooperative effect: Detun- 



ing a single unit destroys the synchronization of the whole 
network. 

The work of Ido Kanter is partially supported by the 
Israel Science Foundation. 
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